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Abstract 

It is shown that response properties of a quantum harmonic oscillator are in essence 
those of a classical oscillator, and that, paradoxical as it may be, these classical 
properties underlie all quantum dynamical properties of the system. The results are 
extended to non-interacting bosonic fields, both neutral and charged. 
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1 Introduction 



The ultimate goal of tliis and forthcoming papers is to extend the phase-space 
techniques developed initially within the quantum-optical paradigm [l|, |2| to 
the relativistic quantum field theory. To this end, we introduce the concept of 
quantum- statistical response formulation of an arbitrary quantum-field system 
and show its equivalence with the conventional Green-function techniques of 
the quantum field theory We prove that Green's functions in the quan- 

tum field theory express, in a somewhat tangled manner, quantum-statistical 
response properties of the quantum system. In short, we show that the ac- 
tual content of the quantum field theory is microscopic signal propagation. 
Response formulation employs the phase-space concepts making it a natural 
"bridging" approach between the quantum field theory and the phase-space 
techniques. We start our analyses with the simple case of a harmonic oscilla- 
tor. Results found for this generic system then serve as leading considerations 
when analysing response properties of more complicated quantum systems. 
In this paper, we show that response formulation may be introduced for any 
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noninteracting bosonic system, whether neutral or charged. In forthcoming 
papers we further extend the results to arbitrary interacting quantum fields, 
both bosonic and fermionic. 



Response is a universal property of any physical system relating all other 
properties of the system to experimentally observable quantities. In particular, 
formulating a system in terms of its response removes the annoying incom- 
patibility between the quantum and classical techniques (q-numbers versus 
c-numbers). A system that is physically classical (i.e, non-quantum in its ob- 
servable properties) becomes in a response formulation also formally classical. 
Response reformulation of a physically classical system recovers a c-number 
description of the same system. This occurs the facts notwithstanding that this 
system has been initially formulated in q-number terms, and that the overall 
approach remains fully quantum. In this paper we show in detail how this 
"recovery of classicality" in quantum mechanics emerges for the simplest ex- 
ample of a harmonic oscillator. We demonstrate that response properties of a 
harmonic oscillator are exactly those of the classical oscillator, and that, para- 
doxical as it may sound, these classical response properties, in fact, underlie 
all quantum dynamical properties of the oscillator. A classical interpretation 
may be given even to the vacuum fluctuations in the harmonic oscillator which 
are commonly believed to be the quintessence of the quantum. 

A word of caution is in place here. The literature on the harmonic oscillator 
is enormous, and majority, or may be even all, of the technical points we 
address in this paper are included in the earlier work, see, e.g., @,0,i. What 
we regard as the actual result is not these formulae as such, but the pattern of 
response fomulation we introduce based on them. The utility of our approach 
will eventually emerge when we extend the treatment to interacting fields in 
forthcoming papers. 



Perhaps the best way of putting our results, however informally, into the right 
historical perspective is, "Schwinger=Kubo+Glauber." Schwinger recog- 
nised the importance of formal c-number currents in the equations for Green's 
functions in the quantum field theory. Kubo \IQ\ formulated his linear response 



theory by adding classical source terms to quantum Hamiltonians. Glauber [11 



was first to understand the relevance of phase-space concepts to macroscopic 
measurements. Our analyses show that Schwinger's c-number currents should 
be kept distinct from Kubo's external sources. This is especially clear within 
the so-called closed-time- loop formalism by Schwinger lij, see also [5|. In this 
formalism, a pair of c-number Schwinger currents are introduced governing 
quantum evolution forwards and backwards in time. We show how one can 
rearrange these into two combinations, the first corresponding to the Kubo 
source, and the second — to the self-radiation of the system in the presence of 
this source. Full information on response properties of the system hidden in 
Schwinger's closed-time-loop formulation is then conveniently reexpressed by 
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combining Kubo's and Glauber's approaches. Detailed discussion of these and 
related points will be given in section HI 



The paper is organised as follows. Quantum-statistical response properties of 
the oscillator are analysed in section[2l Extension of these analyses to arbitrary 
noninteracting bosonic fields is given in section [31 In section 14.11 we introduce 
a response formulation of noninteracting bosons and clarify its relation to 
Glauber's photodetection theory 0, Q Q Q. In section we discuss the 
relation between the response formulation and Kubo's linear response theory 
lo| . In section 14.31 we compare the concept of external c-number current in 
Kubo's and Schwinger's approaches. In appendix [A] we outline Hori's form of 
Wick's theorem 15| in the closed-time-loop formalism. In appendix[B]we prove 
that an alternative way of defining the Kubo current through the Schwinger 
currents discussed in 14.31 is indeed related to the symmetric operator ordering. 



2 Quantum-statistical response of the harmonic oscillator 



2.1 Harmonic oscillator basics 



We consider the quantum harmonic oscillator, with the Hamiltonian in the 
Schrodinger picture 



"2 2 "2 



hujn ( a' a + - ) . 



Here q and p are the standard position and momentum operators 



[q,p] = ih, 



and a^, a are the usual creation and annihilation operators, 



a, a 



h 



muQ 



Po = JhrnuQ. 



(3) 
(4) 



Furthermore, the vacuum (ground) state of the oscillator is defined by 



a|0) = 0. 



(5) 
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Assuming that the Schrodinger and Heisenberg pictures coinside at to = 0, it 
is easy to find the Heisenberg operators related to the Hamiltonian ([1]): 



a{t) = ae-*"°*, a\t) = a^e^'o*, 



m 

Pit) 



V2 
V2 



a{t) - a\t) 



(6) 



where a{t) = e i^/^)Hot^Q(t/h)Hot^ ^^^^^ similarly for q{t) and p{t). Finally, the 
normal ordering places all a^s to the left of all ds, e.g. 



:q{t)q{f): = ^ [a(t)a(t') + a\t)a{t') + a\t')a{t) + a\t)a\t') 



(7) 



2.2 Driven oscillator 



Following Kubo [10], we analyse response of the oscillator by adding a source 
term to the Hamiltonian ([1]) 

ns{t) = nos + j{t)q. (8) 

For simplicity we assume that j{t) =0, t < 0. We consider the source term 
as a perturbation, but do not assume it to be small. In this terminology, T^os 
is a free Hamiltonian and Ti-sit) is a full Hamiltonian, while q{t) and p{t) are 
the interaction-picture operators. The Heisenberg position and momentum 
operators related to Tisit) will be denoted as qj{t) and Pj{t). The presense of 
the source term in ([8]) makes transition to the Heisenberg picture nontrivial. 
The Schrodinger and the Heisenberg representations of an arbitrary operator 
{xs{t) and Xj(t), respectively) are coupled by the unitary evolution operator, 

Xj{t) =U\t)xs{t)l({t). (9) 
The evolution operator obeys the Schrodinger equation, 

ihU{t) = nos{t)U{t), (10) 
while Xj{t) obeys the Heisenberg equation 

ihijit) = \xj{t),nj{t)] +U\t)isit)U{t), (11) 
where iijit) is the full Hamiltonian in the Heisenberg picture, 

nAt)=uKm,smt) = 1^ + '^^^^^ + mqAt)- (12) 
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In particular, 



ihpjit) = \pj{t),nj{t) 



(13) 



Commutators in f|T3l) are easily calculated resulting in the classically-looking 
pair 



m 



(14) 



Substituting Pj(t) = mqj(t) into the second equation results in the inhomoge- 
neous "wave" equation for the oscillator displacement 



m 



m 



while q{t) obeys the homogeneous equation 

'q{t)+uj'om = 0. 
By noting that qjit) = q{t) for t < we find 

q^{t)=q{t) + q,{t), 

where 

q,{t) = J dt'Dn{t-t')j{t'). 



(15) 



(16) 



(17) 



1^ 



Omitting limits of a time integration means that this integration is from plus to 
minus infinity, cf. endnote 19|]. In (|T8|) . is, up to overall sign, the retarded 
Green's function of the c-number equation for the displacement of a classical 
oscillator under the influence of an external force f{t) = —j{t) 



q{t) + ujlq{t) 



fit) 



m 



q{t) = qin{t)- Jdt'D^it-t')f{t'), 



(19) 



where qin{t) is a solution with j = 0. The kernel obeys the typical inho- 
mogeneous "wave" equation for Green's function. 



Z)R(r)+a;2DR(r) 



6{t) 



m 



Explicitly, for the harmonic oscillator ([T]) we have 



sm ujqT 



(20) 



(21) 
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Equation (fT7|l has an important consequence. The difference between qj(t) and 
q{t) is a c-number and does not interfere with operator orderings. The concept 
of normal ordering defined for the interaction-picture operators q{t) and p{t) 
thus equally apphes to the Heisenberg-picture operators qj{t) and Pj(t). 



2.3 Quantum-field-theoretical approach to the driven oscillator 



The goal of our analyses is to find a way of approaching interacting non- 
linear quantum systems. To this end, we treat the oscillator using concepts 
which are normally reserved for much heavier problems, where a solution in 
terms of operators is a rare exception. Namely, we enaploy the well-known 
closed-time-loop Schwinger-Perel-Keldysh formalism js], applicable to any 
nonlinear nonequilibrium quantum-statistical problem (see, e.g., a review by 
Fred Cooper 16|). In this formalism, the oscillator is characterised in terms 
of the double-time- ordered averages of the displacement operators 



{T_q{t,)---q{t^)T+q{t\)---q{t'^)) 



5r^^{t^)---5r^.{tm)5Ti+{t',)---5r^+{t'^) 



r)-=r]+=0 



where we have introduced the generating, or characteristic, functional 



(22) 



= 1 - 



T_ exp 



E 

m+n>l 



i J dtr]^{t)q{t) 



mw. 



T+exp ~ J dtr]^{t)q{t) 
X (T_g(ti) ■ ■ ■ qit^)T+qit[) ■ ■ ■ q{t'^)) , (23) 



)m+n 



with the functional arguments ri±{t) being arbitrary smooth c-number func- 
tions. The symbol (/ dt)^^^ denotes integration over all time arguments, cf. 
also endnote [loj. Time-ordering T+ (T_) puts operators in order of increasing 
(decreasing) time arguments. That is, if ti < t2 < ■ ■ ■ < tm, 



T_g(ti) ■ ■ ■g(t^ 
T+q{ti) ■ ■■q{tr^ 



q{ti) ■ ■ ■q{tm-i)q{tm), 



(24) 



By definition, bosonic operators under the ordering sign commute. The quan- 
tum averaging in fl22|) . fl23|) is defined in the standard way with respect to the 
initial state of the oscillator, described by the density matrix po; 



) = Trpo(---)- 



(25) 
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Fig. 1. C-contour (thin lines) travels first from t = —oo to t = +00 (direct branch) 
and then back to t = —00 (reverse branch). All operators on the direct branch 
are placed to the right of all operators on the reverse branch, while operators on 
the direct (reverse) branch are (!"_) ordered among themselves. This ordering 
rule obviously reproduces the double-time ordering defined by Eq. ([22|) . The three 
contractions defined by Eqs. (I38p - (l40p are shown schematically as thick lines. They 
correspond to three possible ways for an operator pair to be positioned on the C- 
contour: both operators on the direct branch, both operators on the reverse branch, 
and one operator on the direct branch, the other on the reverse branch. 

Functional (!23|) obeys the reality property, 

^*{v.,r]+) = Hvl,V-), (26) 

which is a consequence of 

[T^x{t)---mf = T^S:^it)...y\t'), (27) 

where x{t) ■ ■ ■y{t') are arbitrary operators. This formula expresses the simple 
fact that Hermitian conjugation inverts the order of factors turning a T+- 
ordered product into a T„-ordered one, and vice versa. 

The double time ordering defined by Eq. fl22l) is often understood as operator 
ordering on the so-called C-contour shown in Fig. [H It travels first from t = 
—00 to t = +00 and then back to t = —00 (hence the name "closed-time- 
loop formalism"). All operators on the direct branch are placed to the right 
of all operators on the reverse branch, while operators on the direct (reverse) 
branch are T_|_ (r_) ordered among themselves. This ordering rule obviously 
reproduces the double time ordering defined by Eq. (!22|) — note that "earlier" 
operators go to the right. However, for our purposes it is more convenient to 
keep the T-t-orderings explicit. The interpretation of $(?;_, ?7+) as Schwinger's 
closed-time-loop amplitude jl2| will be discussed in 14.31 

Definitions identical to (!22|) . (!23l) are also given for the quantities related to 
operator qj{t): 



(T_^,(ti)---g,(Ur+g,(t;)---g,(0) 

(-z)'^z" (5™+"<l>(77_,r7+;i) 



r]-=r]+=0 
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where 



T_ exp 



i / dtri_{t)qj{t) 



r+ exp 



i / dtri+{t)qj{t) 



>. (29) 



Equation ( fTTI) then results in a simple formula expressing $(r^_,77+;j) by 
$(r/_,r/+): 



Here, 



and 



?7(t) = -2[77+(t)-r/_(t)], 



'^'cK^; j) = exp J dtr]{t)qj{t), 



(30) 



(31) 



(32) 



where qj (t) given by Eq. (ITSl) is the classical displacement caused by the source 
j{t). The quantity $ci(^;j) is nothing but the characteristic functional for 
products of these displacements at different times: 



Qji^i) ■ ■■qj{tm) 



5™<l>ei(r/;j) 



driiti) ■ ■■5T]{tm) 
Obviously, $ci(^;i) is a purely classical object. 



rj=0 



(33) 



Equation fl30l) is a good example of how manipulating charactestic functionals 
gives access to physics contained in quantum Green's functions. It is also the 
first indication that merely by introducing new functional arguments one can 
make this physics much more transparent. 



2.4 Wick's theorem and vacuum Green's function 



Relation fl30|) expresses two well-known properties of the harmonic oscillator: 
that the classical external source creates a coherent amplitude, and that in- 
trinsic quantum properties and response of a linear system are independent, 
hence the factorisation. To calculate $(?]_, //+), we make use of another case 
of quantum-statistical independence and ensuing factorisation characteristic 
of linear systems. Namely, 

^V-,V+) = ^vUV-,V+)^in{v), (34) 
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where 



$i„(r7) = ^:exp J dtr]{t)q{t): 



(35) 



characterises the initial state of the oscillator, and 



^'vac(^?-,?7+) 







r_ exp 



i / dtri^{t)q{t) 



T+ exp 



t / dt7]+it)q{t) 



0), (36) 



is the functional $(//_, ?7+) for the oscillator in a vacuum state. The functional 
variable ^^(t) in flM|) is the same as appearing in Eq. fl5U]) and defined by flM]) . 



To prove equation fl51|) it suffices to note that it is nothing but a closed form 
of Wick's theorem for the double-time ordering. To clarify this connection 
we look for the explicit expression for $vac('7-, ??+)• Referring the reader for 
details to appendix |Al here we proceed by noting that the harmonic oscillator 
in a vacuum state is a Gaussian system, and, as for any Gaussian system, the 
functional $vac must be an exponent of a quadratic form. 



^vaciV--^V+) = exp < ih J dtdt' 



--r)+it)r]+it')Dp{t-t') 



+ - V-{t)V-it')DUt - t') + V-{t)v+it')Dit - t') 



(37) 



where notice was also taken of Eq. fl2Bl) . It is easy to see that the kernels here 
are the "contractions," i.e., the differences between the time and normally 
ordered products of pairs of operators. Indeed, by using Eq. (l37j) to calculate 
pair averages we have: 



ihDpit - t') 



?7+=r?_=0 

(0 \T+mm\o) = T^mm - ■mm-, m 



- ihD*p{t - 1') 



^^$vac(?7-,?7^ 



5ri^{t)Sr]^{f 



r}+=V-=0 

(0 |T_g(t)g(t')| 0) = T_g(t)g(0 - :g(t)g(t'):, (39) 



ihD{t - t') 



57]4t)6r]4t>) 



r]+=ri-=0 

= {o\mm)\o) = mm)---mmy- m 
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The contractions, which in the quantum field theory are often called free-field 
Green's functions, are expressed by commutators and are therefore c-numbers, 
while their expressions as vacuum averages follow from (0 \:q{t)q{t'):\ 0) = 
(which holds for any normally ordered product). In the Feynman diagram 
techniques [s], Q] -Df emerges as a propagator, while Dp and D occur also as 
parts of the matrix propagator in the Perel-Keldysh techniques jil . By making 
use of Eqs. (Q, for the harmonic oscillator we find 

D{t) = , Dp{t) = e{t)D{t) + e{-t)D{~t). (41) 



We see that equation (IMl) has "Wick's theorem" written all over it. Its LHS 
is a characteristic functional of the double-time-ordered operator averages. Its 
RHS contains the contractions and normally-ordered operator averages. The 
three contractions defined by Eqs. (138|) - (H0|) correspond to the three ways for 
an operator pair to be placed on the C-contour: both operators on the direct 
branch, both operators on the reverse branch, and one operator on the direct 
branch, the other on the reverse one, cf. Fig. [H On expanding both sides 
of (1341) into functional Taylor series and comparing the coefficients, we find 
expressions for the double time-ordered operator averages in terms of sums of 
the normal averages with contractions, just as stipulated by Wick's theorem. 
The remaining tedious task of verifying that all terms in these sums indeed 
occur with coefficient 1 may be replaced by a direct algebraic proof of ( !34l) 



based on Hori's form of Wick's theorem [l5l]. For details of this proof we refer 
the reader to appendix |Al 



2.5 Are vacuum fluctuations quantum? 

2.5.1 Algebra of Green's functions 

By combining (I5D]) with fl5^ we find 

^{r]_,r]+]]) = $vac(r/-,??+)$in(??)$ci(r/;j)- (42) 

Quantum properties of the driven oscillator thus comprise three independent 
contributions: the vacuum fiuctuations (characterised by $vac), the initial state 
(characterised by $in), and the external source (characterised by $ci)- Of the 
three contributions, that of the external source is fully classical. Contribution 
of the intial state is quantum or classical depending on whether this state is 
quantum or classical. As to the vacuum fiuctuations, these are commonly re- 
garded as the quintessence of the quantum. However, contrary to the common 
wisdom, we are going to show that $vac allows for a fully classical interpreta- 
tion in terms of the response of the classical oscillator. 
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As a first step to demonstrate this, we express Dpit) and D{t) by D^{t). It is 
easy to verify tliat tlie kernel Dji{t) in tlie quantum treatment may be defined 
as the retarded combination 



D^{t) = Dpit) - D{-t) = 9{t) [D{t) - D{-t) 
From this equation we find 



(43) 



Dnit) - D^i-t) = e{t) [D{t) - D{-t)\ - e{-t) [D{-t) - D{t) 

= D{t)-D{-t). (44) 



Physics behind Eqs. f HSl) and fj44l) is discussed in section 1421 Here we proceed 
formally by noting that D{t) may be expressed by D^lt) if retaining only 
half the frequency spectrum of Eq. 0441) . Indeed, in the terminology of quan- 
tum optics, D{t) is purely frequency-positive, D{t) oc e~*'^°*, so that D{—t) 
is purely frequency-negative, D{—t) cx e*'^°*. In general, the frequency-positive 
and negative parts g^'^^t) of some function g{t) are defined by dropping the 
corresponding half (negative or positive) of its Fourier-spectrum 



g'^Kt) 



^e-^-%±co)g^, g^ 



dte"^'g{t). 



Obviously, 



g^^\t)+g^-\t) = git). 
This relation, however, implies that g^ is smooth at u = 0. 
Retaining only the frequency-positive part of (jSj) yields 



(45) 



(46) 



D{t)=Di:\t)-D\^>{-t), 



(47) 



where D^^\t) and \t) denote the frequency-positive and negative parts of 
D^{t), respectively. In obtaining Eq. P7|) we have used the fact that invert- 
ing the time argument swaps the frequency-positive and negative parts of a 
function. Finding Dp is now a matter of straightforward algebra. Again using 
T3l), this time to express Dp, we have 



Dp{t) = Dp,{t) + D{-t) = Dn{t) + D'^\~t) - D'^>{t) 



Di^\t)+D^^>{-t), (48) 



where use was also made of Eq. (146!) . 



2.5.2 Response substitution 

In the common parlance, Eq. (1371) embodies vacuum fluctuations in the oscil- 
lator. However, in view of relations (H71) and (HSI) . how quantum are these? 
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We note that is real; for a real function, complex conjugation swaps its 
frequency-positive and negative parts, 

DW*(t) = DW(t). (49) 

Complex-conjugating fH8l) then allows us to express Dp{t) occuring in fl37|) as 

DUt) = D^n\t) + D^^\-t). (50) 

Consider now the first term in Eq. (1371) . By making use of Eq. (HHI) we have 

1 



ih J dtdt' - -ri+{t)ri+{t')DF{t - t') 

= -| / dtdt'ri+{t)r,^{t') [4+)(t - t') + 4+)(t' - t) 

= -ih j dtdt'r]+{t)D^^\t - t')r]+{t'). (51) 

We now employ the following connection between the frequency-positive and 
negative parts and the convolution: 



dt'h^^\t - t')g{t') = / dt'h{t - t')g^^\t'). 



(52) 



This is proved by noting that the Fourier-transform of either expression is 
9{±u!)g^h^. From Eqs. flSTj) and (l52l) we find 

1 



ih J dtdt' - -r]+{t)r]^{t')DF{t - t') 



ih I dtdt'r]+{t)DKit - t')r]^+\t'). (53) 



Similarly, by making use of Eqs. fH7|) . flSUjl and (15^ . we transform the two 
remaining terms in Eq. (1571) as follows: 



+ 2^-(^)^-(0^F(i-0 
= Y J dtdt'r]_{t)ri_{t') 



d^\t-t') + D'^\t'-t) 



(-)/ 



ih J dtdt'r]_{t)Di^\t - t')r]_{t') = ih J dtdt'r]_{t)D^{t - t')ri^s\t'), (54) 



ih dtdt' 



+ r^_{t)v^(t')D{t-t') 
= ihj dtdt'r]4t)r]+{t') D^^^t - t') - D\^\t' - t) 
= ihj dtdt'7]_{t)D^^\t - t')r]+{t') -ihj dtdt'r]+{t)D^^\t - t')7]^{t') 
ih j dtdt'ri_{t)D^{t - t')ri^_^^\t') -ihj dtdt'ri+{t)D^{t - t')ri^_r\t'). (55) 
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Combining Eqs. (153|) . (IMl) and (l55l) yields 

$vac(^?-,?7+) 



exp |— i/i y (itfit' 



?7+(t)-?7_(t) 



I5R(t - 1') 



■ (56) 



On introducing a new set of functional variables, 
we can write flS^ as 



(57) 



vac (??-,??+) = exp 



dtr]{t) J dt'D^it - t')a{t' 



(5J 



Note that the functional variable rj{t) of Eqs. (!30!) and flMl) has occurred once 
again. 



Equation (1^ . we remind the reader, expresses vacuum fluctuations of the 
quantum harmonic oscillator, and yet it looks astonishingly classical. In the 
variables r/, a the functional $vac has turned simply into a generating functional 
for products of a c-number field emitted by c-number source a: 



vac (??-,??+) = ^c\{ri\(T), 

where $ci was defined by Eq. 



(59) 



Of crucial importance is that fl57p is a genuine change of functional variables 
in the characteristic functional. We show this by solving for the inverse sub- 
stitution. On taking the frequency-positive and negative parts of Eqs. (1571) we 
have 



zry(-) (t) = r/t ^ {t) - r]^S^ (t) , a^-\t) = hr]^'^ (t) . 
Solving for the frequency-positive and negative parts of r]±(t) yields, 

vi^\t) = iaW(t), vi-\t) = zv^-\t) + ia(-)(t), 

n n 
By using Eq. fj46|) . the inverse of fl57j) is found to be 



(60) 



(61) 



(62) 
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proving that ?7+(t), r7_(t) ^ 'r]{t),a{t) is indeed a genuine functional substitu- 
tion. 



Substitutions (!57|) . (!62|) are inherently complex. If we choose ?^±(t) to be real, 
then 77(t),cr(t) will be complex, and vice versa, but we cannot make all four 
functional arguments real simultaneously. It makes obvious physical sense to 
keep a{t) real, while for the other three it is of little consequence. A convenient 
formal choice is = rj"^, this making 77, a real and also turning $(?7-, into 
a real quantity, cf. Eq. (I23l) . 



3 Response of free bosonic fields 

3.1 Neutral bosons 
3.1.1 Definitions 

In this section we extend the treatment of the harmonic oscillator to gen- 
eral noninteracting bosonic fields. We restrict our attention to the key point: 
derivation of field versions of the response substitution (!57|) . transforming the 
quadratic form of contractions (1371) into the classical expression (1581) . Here, we 
consider neutral bosonic fields. Charged fields will be the subject of the next 
section. 

To be specific, we consider a Hermitian 4- vector field described by the free-field 
operator Q^{r, t). This choice suggests quantum electrodynamics. However, all 
our manipulations apply to the time variable and do not touch field "labels." 
In particular, we do not make any use of the fact that // is a 4-vector index. 
Below, /i may be regarded as an arbitrary index, e.g., a mode index running 
over a given set of modes. This allows all formulae to be easily adjusted to the 
case of an arbitrary bosonic field. 

Properties of the free-field bosonic operators are considered in any text book on 
the quantum field theory, see, e.g., [3|. The standard quantisation approach 
to the bosons is based on replacing the c-number amplitudes in the mode 
expansion of the classical field by q-numbers: 



Here, Q^^{r) are mode functions, uok are mode frequencies, and a\,ak are the 
standard mode creation and annihilation operators. 



Q,{r,t) = h}/^Y1 a,e— ^*gJ(r) + aIe-^*Qj*(r) . 



(63) 



1 - 



(64) 
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We do not consider normalisation conditions for modes, allowing us to include 
all frequency-dependent coefficients in (1641) into the mode functions. The over- 
all factor h}/"^ in Eq. fl63p allows for a convenient classical limit where h^^'^ak 
is replaced by the classical amplitude. It also makes the mode functions es- 
sentially classical quantities (e.g., defined by classical boundary conditions). 

With generalisation to interacting fields in mind, we wish to eliminate explicit 
mode expansions from considerations. To this end, we now summarise the 
properties of the free-field operators essential for our purposes. The free-field 
operator is naturally split into the sum of its frequency-positive and negative 
parts. 



g^(r,t) = g(+)(r,t) + g(-)(r,t). 



it 



(65) 



The frequency-positive and negative parts always apply to the time variable. 
Obviously, 



gW(r,t) = ;iV2^a,e— ^*Qj(r), 

k 



(66) 



Their commutator is a c-number, 



Q^;\r,t),Ql7\r',t': 



while two Q^'^^s and two •'s commute: 



0, 



(68) 



cf. Eq. ^ for the oscillator. As seen from Eq. (1^ . D^^i{r, r', t) is frequency- 
positive in respect of its time argument. 



Di-J{r,r',t)=0. 



(69) 



The normal ordering is defined by placing all ^s on the left of all Q^~^^s: 

= Ql-\r',t')Qi^Hr,t), (70) 



and similarly for a larger number of factors. 
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3.1.2 Algebra of free-field Green's functions 

We introduce the three contractions in full analogy with the oscillator case: 

ihDF^^,{r,r' ,t - t') = T+Q^{r,t)Q^,{r\t') - :Q^{r,t)Q,,{r',t'):, 
-ihD;^^, (r, r', t-t')= T_Q,ir, t)Q^,{r', t') - :Q,{r, t)Q^.{r\ t% (71) 
ihD,,^,{r, r', t-t') = Q^ir, t)Q^\r\ t') - :Q^{r, t)Q^\r\ t'):. 

Obviously D^^r{r, r', t) here is the same as in Eq. fl^ . It is also easy to check 
that 

DF^,^,{r,r',t) = 9{t)D^^,{r,r',t) + e{-t)D^,^{r',r,-t). (72) 
Similarly to Eq. (143|) we define the retarded combination of contractions 

DR^,^,'{r,r',t) = DFf,f,>{r,r',t) - D ^{r' , r , -t) 

= 9{t) [d^,, (r, r', t) - D^,^{r\ r, -t)] . (73) 

The physical meaning of this quantity will be discussed in detail in the next 
section. By making use of fl73|) we recover the field counterpart of Eq. 



D^^, (r, r', t) - D^,^{r', r, -t) = D^^^. (r, r', t) - Dn^,^{r\ r, -t). (74) 

We now recall that D^^i (r, r', t) is purely frequency-positive and hence D^i^{r' , r. 
is purely frequency-negative. By taking the frequency-positive part of ([71]) we 
find 

D,,,{r, t\ t) = Dg,)^,(r, r', t) - ^^^(r', r, -t). (75) 



By using equation (1731) to express Dp we then obtain 

D^,,ir, r', t) = ,(r, r' , t) + Dit)^{r', r, -t). (76) 

Equations (1751) and (|7B]) generahse (H7|l and to the case of a neutral 
bosonic field. We see that the mode expansion enters only through equation 
fl671) for the D-function. If we assume the D-function, or which is the same, 
the field commutator, known, the concept of mode is eliminated altogether. 



3.1.3 Response substitution for free bosonic fields 

The key result on which the response reformulation of the oscillator is based 
is the transformation of the quadratic form of contractions into the classical- 
emission form, cf. Eqs. (lS7j) and (ISSjl . For the field case, this means representing 
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the quadratic form of the field contractions (cf. Eq. (137|) ) 

In $vac (??+,??-) 

ifiJZj dtdt'd^rSr' 

+ \ri'- (r, ml (r', t')D*^^^, (r , r', t - t') + (r^, (r', t') , r', t - t') 

(77) 

as 

In $vac(^+,??-) 

= ^ 1 dtdt'd^rd^r'r]''{r, t)Dn,,'{r, r', t - t')a^'\r' , t'), (78) 
mm' 

where D-^fj,^i{r,r\t — t') is the retarded combination of contractions defined 



by (1731) . In omitting integration limits in (1771) and (1781) we follow endnote 19 
The relation between r]± and 77, a is found by substituting Eqs. (1751) and (1761) 
into (I77j) and manipulating the emerging expression in detailed analogy with 
Eqs. ( l5T]) - (l55l) . As a result we indeed arrive at Eq. (I78|) with 



?7^(r,t) = -z[<(r,t)-?7^(r,t)J, 
a^(r,t) = ;i[r/';(+)(r,t) + ?7^(")(r,t)]. 

This is a field counterpart of the substitution (1571) . It is now clear that any 
relation for the oscillator may be generalised to the field case by simply adding 
the field "labels" to the time argument, f — >■ /i, r, t. It may be easily checked 
that this indeed generalises all relations derived in the previous section to 
neutral bosonic fields. We also note that substitutions (I79|) bear no trace of 
the mode expansion (|63|) which we used as a starting point. 



3.2 Charged bosons 



3.2.1 Definitions 

Referring the reader for details to any textbook on the quantum field theory, 
here we only list the key properties of the free-field operators. A charged 
bosonic field is described by a conjugated pair of field operators, 

F^(r,t) = iW(r,t) + S(-)(r,t), 
Fl{r^) = M-\r,t) + B^;\r,t). 

We use distinct notation for the frequency-positive and negative parts of 
the field operator to emphasise their distinct physical nature: F^+^(r,t) = 
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A^^\r,t) annihilates particles while F^~)(r,t) = B^~\r,t) creates antiparti- 
cles. Similarly, F^*^+''(r, t) = Bj^^r, t) annihilates antiparticles while F^^^\r, t) 
Ajj^\r,t) creates particles. The operators A^^\ A^~\ B^~^\ and B^~^ pairwise 
commute except the two pairs: 

AM(r,t),Al-\r',t')] 
'_Bl^\r,t),Bl;\r',t') 

Function is frequency-positive while is frequency-negative, 
Df^t\r,r',t) = D;;'^:\r,r',t)=0. 



hDX{r',r,t'-t). 



(82) 



3.2.2 Algebra of Green's functions 

From now on, for brevity we drop the field labels; similarly to the neutral case 
they can always be restored by replacing t —>■ fj.,r,t. Unlike in the neutral case 
we have now four contractions, 

T+F{t)F\t') - :F{t)F\t') 
T^F{t)F\t') - :F{t)F\t') 

F{t)F\t') - :F{t)F\t') 

F\t')F{t) - :F{t)F\t') 

where the normal operator ordering is defined in the usual way as 

= :A^-\t')A^+\t): = A^-\t')A^+\t), (84) 

and similarly for the Bs. Hermitian conjugation of c- number kernels is defined 
in the standard way, 

DUt) = D*{-t). (85) 

This formula implies that we regard Dp as a c-number operator kernel Dp {t — 
t'), and similarly for other kernels. The second line in (!83|) is Hermitian con- 
jugation of the first one. The other two expressions are independent; they are 
Hermitian so that D"^ and D^ are anti-Hermitian: 

D^^t) = -D^{t), D^'^t) = -D^'it). (86) 

The time and normal orderings apply formally to commuting pairs as well but 
in fact do nothing, e.g. 

T^F{t)F{t') = :F{t)Fity. = F(t)F(t'), (87) 
etc. The corresponding contractions are zero. 



ihDp{t - t'), 
-ihDlit-t') 
ihD^{t-t'), 
ihD^{t-t'), 



^3) 
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It should not be overlooked that Dp{t) is now not symmetric with respect to 
time reversal. From (IHT]) and (!83|) . 



D,it) = eit)D\t)+9i-t)D^it), 

DUt) = e{-t)D''^{t) + e{t)D''^{t), ^ ^ 

where the second line is conjugation of the first, cf. flS^ . It may therefore 
seem that one can construct two independent retarded combinations: one by 
subtracting from Dp and the other — by subtracting D"^^ from Dp. These 
turn out to be two ways of expressing the same function 

DR(t) = DfW - D'^it) = e{t) [D^it) - D\t)\ , 
D^{t) = Dlit) - D^\t) = e{t) [D^t(^) _ D^^t)] . 

Equivalence of these two definitions is clearly seen from (1861) . We note that 
Dr is now complex so that one can indeed say that we have two real retarded 
combinations, but this does not appear to have any physical implications. 
Having a complex retarded Green's function is a natural property of a charged 
field which is a quantised equivalent of a complex c-number field. 



Combining the definitions of we find the "charged" equivalent of 

D^{t)-Di{t) = D\t)-D''{t). (90) 

By using (182!) we find 

D^{t) = Dt\t)-D^^\t), 
while by reusing the definitions of to express D-p and we have 

In obtaining the second equation here we used that 

Dr\t) = D^^^\t). (93) 

This holds for any kernel; it also makes it evident that the two equations in 
(p2l) are indeed conjugate. 

3.2.3 The response substitution 

By making use of Eqs. (!9T|) and (!92l) . one can identify the response substitution 
for charged fields. In full notation the "charged" version of Eq. (1781) is written 
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as (cf. endnote |l9l ]) 



ihJ2 J dtdt'd'rd'^r'\-r]^l{r,t)r]i{r',t')Dp^^ir,r',t-t') 
mm' 

J2 J dtdt'd'rd'r' y{r, t)Dj,,^,{r, r', t-t')a^'ir', t')+v^{r, t)D*^^^,{r, r', t-t')a'^' {r' , t') 
mm' 



where 



r]^{r, t) = —i r]^{r, t) — r]^{r, t) 



fj'tir,t)-f]^lir,t) 



a^(r, t) = h[f]'}^^ (r, t) + f]-^'\r, t) . 
Equations (^5^ are simply the substitution (179]) repeated twice. 



(95) 



(94) 



4 Results and discussion 



4-1 Response reformulation of bosons 



Having proven that the results for the harmonic oscillator are trivially gen- 
eralised to arbitrary noninteracting bosonic fields, for simplicity we return in 
the ensuing discussion to the oscillator case. From equations ( H2|) and (!59|) we 
have 



(96) 



We see that the variable a is redundant: all quantum dynamical and statistical 
properties of the oscillator emerge as its self-radiation and response properties 
through the functional 



:exp I dtr]{t) g(t) + J dt' Dn{t - t')j{t') (97) 



Describing the system in terms of the functional $R(?7;j) will be called a 
response formulation. This functional can be re-expressed in two equivalent 
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forms. Firstly, by noting that the quantity in square brackets equals q{t) + 
qj{t) = qj{t), cf. Eqs. (fT7|) and (fTSi) . we can write 



^r{v;j) = (:exp J dtr]{t)q,{t):'^ . (9J 



That is, ^ulrj] j) is a characteristic functional of the normal averages of oper- 
ator qj{t). Equation flUSl) thus establishes a connection between the response 
formulation and Glauber's macroscopic photodetection theory 0, 11, 13, 14 



Suitably generalised, this relation will play a fundamental role in the analyses 
of interacting systems. Secondly, introducing the P-distribution P{a) charac- 
teristic of the initial state, 

(:at"a™:) = [ ^P{a) a^'^a"' = (99) 



we can rewrite fl971) as 



$R(r7; j) = exp J dtr]{t) q^t) + J dt'Dnit - t')]{t') 



(100) 



where 



gi„(t) = -^(ae-^'^°* + aV-o*). (loi) 

We use the bar as a notation for the averaging over the initial P-distribution 
to emphasise the analogy between the response formulation of the quantum 
oscillator and a c-number description of the classical oscillator. Were we to give 
this symbol back its customary meaning of a classical statistical averaging, 
f llOlj) would become a characteristic functional of multi-time momenta of the 
c-number oscillator coordinate q{t): 



q{t) = qUt) + I dt'Dnit - t')j{t') 



q{ti) ■ ■■q{tm) 



5'"$r(^;j) 



(102) 



ri=0 



Except for the possible nonpositivity of P{a), Eqs. (11 001) and fll02p appear 
fully — and astonishingly — classical. 

Using ri{t), a{t) as arguments in the charactristic functional brings out perfect 
physical sense. Function a{t) has turned out to be combined with the exter- 
nal c-number driving force, a property characterising the input of a system; 
we shall in future see that this property is quite general and also holds in 
the case of nonlinear interacting fields. Setting a{t) to zero turns $(?;_, r^+;j) 
into a characteric functional for the normally ordered averages; again, this is 
quite natural a property for the output of a system, recalling that the nor- 
mally ordered operator averages emerge in the macroscopic detection theory 
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of Glauber ll|]. We see also that, whether the initial state of the quantum os- 
cillator is quantum or classical, its quantum dynamics, i.e., everything related 
to the response to the driving force as well as vacuum fluctuations, emerges 
in fully classical form. 



4-2 Kubo's linear response as a quantisation condition 



Equation (p8l) makes evident the connection between the response formulation 
and Glauber's photodetection theory. Our immediate goal is to clarify another 
important connection, namely, between the response formulation and Kubo's 



linear response theory [lOj. The way was introduced into the quantum 
treatment in section [2.5.11 gives a misleading impression of pure luck: we just 
happened to notice that a particular combination of the quantum Green's 
functions coincides with the retaded Green's function of the classical oscillator. 
In fact the situation is far more physical. By using the primary definitions of 
Dp and D as contractions, cf. Eqs. fl55]) - fHUl) . we rewrite Eq. fH51) as 

Dnit) = -I [T+mm - mm] = -^oit - n nm, m)]) ■ (103) 

The commutator on the RHS here is a c-number, so setting it formally under 
the quantum averaging does not change anything. We introduced the averaging 
to make it obvious that the RHS of fllOSp is in fact Kubo's formula for the 
linear response function of the quantum oscillator with respect to the external 
current in Hamiltonian ([8]). Thus the actual physical message of equation (H3!) 
is that the linear responses of the classical and quantum oscillators coincide. 
In other words, if we restrict the measurement to the linear response, the 
quantum and the classical oscillators become experimentally indistinguishable. 

Statements like this one are at the core of the very concept of quantisation. 
What it says, in essence, is, "the change from the c-number to the q-number 
description preserves certain physical aspects of the situation." The standard 
quantisation approach for the oscillator says, "retain the equations of mo- 
tion for the coordinate and momentum while making these quantities non- 
commuting." Unlike equations of motion, response is an observable property 
so specifying it rather than the equations looks as much more physical an ap- 
proach. We now show that a consistent quantisation approach to the oscillator 
can indeed be devised based on Eq. (11031) . To start with, we show that this 
equation allows one to express the commutator by Dr. Firstly, we rewrite it 
as 

e{t - t') [q{t),Q{t')] = ihDnit - t'). (104) 
Swapping here t ^ t' and then using the antisymmetry of the commutator 
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yields: 



9{t'-t)[m,m\ = -oit'-t)[m,m 

Subtracting fllOSp from (11041) we find 



ihD^{t' - t). 



ih 



Dn{t - t') - Dnit' - t) 



(105) 



(106) 



We see that, for the free field, the information contained in the two-time 
commutator is the linear response of the field. Paradoxical as it may sound, the 
noncommutivity of field operators, this quintessence of the quantum, for the 
oscillator is simply a technical way of describing response. This statement is 
fully supported by the formal treatment and must therefore be taken seriously. 
We shall see in the future that the technical equivalence of noncommutivity 
and response is true also for interacting nonlinear systems. 

Full set of observables for the oscillator includes momentum p{t). Commuta- 
tors involving p{t) are found on supplementing (11061) by 



p{t) = mq{t). 
In particular, by making use of Eq. (PT]) . 

d r. 



(107) 



ihm— 
at 



Dnit-t') - D^it' -t) = ih cos uo{t-t'). (108) 



The canonical quantisation condition follows: 



ih. 



(109) 



We note that the definition of momentum is quite general while all physical 
particulars about the oscillator enter through the response. 



4.3 C-number sources in quantum mechanics: Kubo versus Schwinger 



4-3.1 The Kubo and Schwinger currents 

It is instructive to clarify the relation between the concepts of external c- 
number current in Schwinger's and Kubo's approaches. To this end, consider 
the evolution operator in the interaction picture. It obeys the equation 



ihU{t)=j{t)q{tpit)- 



(110) 



solved by 



Uit) = T+ exp 



h 



dt'j{t')q{t') 



:iiii 
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Following Schwinger 12], we distinguish evolution forward in time from evo- 
lution back in time. The former is determined by current j = j+{t), the latter 
is commanded by an independent current j = j-{t). Formally, this means 
considering a conjugate pair of evolution operators: 



= r+ exp 



dt'uimt') 



(112) 



fl J —CO 

propagates the system forward in time, while the conjugate evolution operator 

UUt) = T.exp \l 
in 



dt'j-{t')m 



(113) 



propagates it back in time. Functional $(?]_, ?7+) can then be written as a 
Schwinger-style amplitude, 



$(r7_,r7+) = (5l5+), 



(114) 



where 5+ and <S_ are the forward-in-time and back-in-time S-matrices, 



5+ = W+(oo) = r+ exp 
Sl = Wl(oo) = T_exp 



h 



h 



dt'Mmt') 



dt'Mmt'] 



(115) 



and 



J±{t)=hr]±{t). 



Transformations fl57ll then read 



hr]{t) 



u{t)-j-{t)], j = jW(t)+ji-)(t). 



(116) 



(117) 



We see that the dependence of the Schwinger amplitude on the j± pair con- 
tains full information on the response properties of the system, but in a fairly 
coded and tangled manner. If = j_, the Schwinger amplitude is simply 
unity. We have to consider j+ and j_ as independent quantities; in suffices to 
assume that j+ is complex and j_ = j^. There is no possibility to have j, j+ 
and j_ real simultaneously. Unlike Kubo's current, Schwinger's currents are 
formal quantities. Physically, it is natural to consider the 'r],j variable pair 
as primary, and the j± pair — as secondary. In other words, we regard the re- 
sponse viewpoint as a primary physical viewpoint. Description of the system 
in terms of its dependence on the Schwinger currents is from this perspective 
secondary. 
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4-3.2 Response formulation and the operator ordering 

One may wonder what if we used some other way of defining the Kubo current 
by the Schwinger currents. For instance, let us define, 

hvit) = -I [3+{t) - 3-{t)] , Jw(t) = \ [3+{t) + 3-{t)] . (118) 

Does functional $ (^^, have any meaning in variables r}.,jY^l Due to ad 
hoc nature of flllSp . it is even unclear if can be interpreted dynamically in 
terms of Hamiltonian ([H]) with j jw ^ priori., calling jvv a Kubo current 
is in no way justified, but this interpretation of indeed holds. Namely, in 
appendix [B] the following relation is derived: 

$(7;_,r/+) = $S'(r/, Jw) = (W^exp j dtr^{t)q,^{t)^ , (119) 

where W ■ ■ ■ denotes the symmetric (Weyl) ordering [2] of the creation and 
annihilation operators. This way, substitution flllSp introduces an alternative 
response viewpoint based on Kubo's Hamiltonian and Weyl's operator order- 
ing. It can also be shown that yet another way of defining Kubo's current, 

jAit)=ji-\t)+£^\t), (120) 

is equivalent to combining Kubo's Hamiltonian with antinormal operator or- 
dering, 

HV-,V+) = ^^kiv.JA) = (^exp I dtri{t)qj^{t)^ . (121) 

It is then only natural to come up with a conjecture that specifying the way 
the Kubo current is expressed by the Schwinger currents equals specifying an 
operator ordering used to define the output of the system. Physically mean- 
ingful discussion of this conjecture must be postponed till general properties 
of response reformulation are understood for interacting systems which is a 
subject of our forthcoming papers. Here we restrict ourselves to two remarks. 
Firstly, we note that this result is certainly counter-intuitive. Contrary to 
expectations, changing the way Kubo's current is expressed by Schwinger's 
currents affects not the fact that Kubo's current may be interpreted as such 
but rather the way we describe the system in phase-space terms. Redefining 
Kubo's current affects not the input but the output of the system. Secondly, 
it is the characteristic property of the normal ordering to map the quantum 
oscillator identically on the classical oscillator. All other orderings introduce 
various "vacuum noises" resulting in apparent dissimilarity of the quantum 
and classical viewpoints. 
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5 Conclusion 



We have introduced a quantum- statistical response formulation of noninter- 
acting bosons relating all quantum dynamical properties of the system to 
its response properties. This formulation is, by construction, equivalent with 
the conventional description of the system in terms of Green's functions in 
Schwinger's closed-time-loop formalism. It is also shown to be naturally re- 
lated to Glauber's photodetection theory and Kubo's linear response theory. 
Its most important property is that a physically classical system is in response 
formulation also formally classical, i.e., described by c-mimbers and probabil- 
ities. In our forthcoming papers we shall generalise the response formulation 
to interacting quantum fields, bosonic as well as fermionic. 
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A Wick's theorem 



In its initial form. Wick's theorem expresses time-ordered products of free- 
field operators by their normal products. It can easily be verified that proof 
of Wick's theorem (see, e.g., (3|) is based solely on the fact that the time axis 
is linearly ordered. It can therefore be extended to field operators defined on 
any linearly ordered set. In particular j^, it holds for the double-time-ordering 
which, as we have already mentioned, can alternatively be regarded as ordering 
on the C-contour shown in Fig. [TJ 



In the case of harmonic oscillator. Wick's theorem for the double-time ordering 
[5I may conveniently be written as a closed relation. 



T_ exp 



i / dtri_{t)q{t) 



T+ exp 



I J dt7]4t)q{t) 
$vac(^-,??+) :exp j dtr]{t)q{t):, (A.l) 



where $vac(^-, "'?+) is given by Eq. (!37|) . Equation (!34l) then follows simply by 
averaging (]A.1|) . We start from proving relation (]A.1|) for the T+-ordering, the 
case of double-time ordering being then a matter of straightforward general- 
isation. Wick's theorem for the T+-ordering stipulates that "a time-ordered 
product of free-field operators equals a sum of normal products resulting from 
all possible ways of replacing pairs of operators by their contractions," the lat- 
ter being defined by Eq. fl38l) . As was first noticed by Hori 15|, the algebraic 



structure corresponding to the "sum of operator products with all possible 
contractions" is produced by an exponential differential operator applied to a 
product of c-number functions q(ti) ■ ■ ■ q{tm)- Hori's observation allows one to 
write Wick's theorem as 



T+q{ti) ■ ■ ■ q{tr 



exp(A) q{ti) ■ --qitr. 



(A.2) 
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where 

A = '- [ dtdt'Dpit - t')—4---. (A.3) 

To prove this, we note that applying A results in correct patterns with one 
contraction: 

Ag(t) = 0, 
/^q{t{)q{t2)=ihDF{h-t2), 
Aq{ti)q{t2)q{h) = ihDpih ~ ^2)^(^3) + (A.4) 
ihDF(ti - h)q{t2) + ihDF{t2 - h)q{t{), 



All terms enter with coefficient 1 as required. Furthermore, produces terms 
with two contractions, each with coefficient 2; this reflects the two ways a pair 
of given contractions can be added sequentially. The same reasoning shows 
that A" produces terms of n contractions with coefficients n\. The correct 
combination of terms with all possible contractions, each entering with coeffi- 
cient 1, including the term without contractions, is therefore indeed produced 
by the sum 

00 -1 

l + E3A" = exp(A). (A.5) 

n=l "'■ 



By making use of flA.2p we have 



r+ exp 



-I / dt7]+{t)q{t) 



: < exp(A) exp 



I / dtri+it)q{t) 



:. (A.6) 



The quantity in curly brackets is easily calculated using that, for any pair of 
functional F{s) and G{q), we have 



q'=0 



q=0 



(A.7) 



This can be proven directly by expanding the functionals into their functional 
Taylor series. By making use of these relations, the calculation reduces to the 
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application of a functional shift operator: 



exp(A) exp 



-t I dtr]4t)q{t) 

ih 

X exp 



exp I ^ y dtrj+{t) 
dtdt'Dpit - t')s{t)s{t') 



Qit) 



s=0 



Ss{t) 



— exp 



i / dtrj+it)qit) 



exp 



^ J dtdt'Dp{t-t')r^+{t)'q+{t') 



(A.8) 



As a result we have 



r+ exp 
= exp 



-I J dtr]+{t)q{t) 
ih 



J dtdt'Dpit -t')r]+{t)ri+{t') :exp -i J dtr]+{t)q{t) 



:, (A.9) 



which is indeed relation ( lA.ll) written for the case of r+-ordering, i.e., with 
r/_ = 0. 



Generalisation of this argument to the double time ordering reduces to apply- 
ing Wick's theorem to the C-contour (Fig. [1]) in place of the time axis, cf., 
e.g., Ref. [sl. Hori's form of Wick's theorem for the double time ordering reads 
(cf. (isl) 



T_q{h)---q{tm)T+q{t[)---q{Q 
exp(Ac) g-(ti) 



■q.{tm)q+{t[)---q+{t'J (A.IO) 



where 



Ac = ih 



J dtdt' 



1 , 52 

2^'^' ~ ''hq4t)Sq^it') 

2^*^^^ ~ *'hq4t)Sq4t') ~ 



D{t - t'] 



6q4t)6q4t') 



(A.ll) 



This can verified in detailed analogy to how Eq. (]A.2I) was proven, cf. Eqs. 
( ]A.4p and (lA.Sp . The rest of the calculation leading from (\A.6\i to (]A.9|) gen- 
eralises trivially resulting in Eq. (lA.ip . 
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B The symmetric mapping 



Here we give a proof of equation flllOp . Regarded as functional variables, 
currents j and jw are coupled by the substitution 



By making use of Eqs. (1961) and (1971) we then find 



(B.l) 



(B.2) 



where 



:exp J dtr]{t)q{t): 

X exp |y I dtdt'r]{t)DR{t - t') r]^+\t') - r]^-\t'] 



(B.3) 



It is a well-known fact (see, e.g., 17|) that characterisrtic functionals cor- 
responding to different operator orderings differ in Gaussian factors. This is 
exactly the case for ^^{r]) and $in(?7), so that ^"^{r]) must be a characteristic 
functional corresponding to some ordering other than the normal ordering. Be- 
low we prove that ^^{r]) is a symmetrically ordered characteristic functional 
of the initial state, 



(B.4) 



Equation (11191) then follows by combining Eqs. flB.2p and flB.4p and recalling 
Eq. (dZD. 



We start the proof of flB.4p from simplifying the Gaussian factor in Eq. flB.3p . 
Equation (152!) allows us to shift the frequency-positive and negative parts on 

J dtdt'^{t)D^{t - t') [r;W(t') - 7^^-\t') 

' dtdt'riit) \D''^\t - t') - d^ \t' - t)] r]{t'). (B.5) 



We have also changed \t — t') D^^ \t' — t); it does not affect the result 
because of the symmetrisation imposed by the quadratic form. Using Eq. (l47l) 
we obtain 



^^iv) = ^:exp J dtvit)qit)-)(^W y / dtdt'r]{t)D{t - t')r]{t') 



(B.6) 
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To prove Eq. (1B.4P it suffices to show that quadratic terms in (1B.4P and (1B.6P 
coincide: 



- J dtdt'r]{t)r]{t') [ + ifiD{t - t') 

1 



- J dtdt'rjitUt') iWrnm) ■ (B.7) 

A nearly identical relation follows directly from the definition of D{t) as con- 
traction by Eq. (HUj) : 



dtdt'r]{t)T]{t') {:q{t)q{t'):) + ihD{t - t' 

1 



dtdt'r]{t)r]{t') (g(t)g(t')) • (B-^ 



In view of the symmetrisation of kernels imposed by the quadratic form we 
are left to prove that 



1 



qmit') + q{t')m] =Wq{t)qit')- 



(B.9) 



This is easily checked by making use of Eq. ([6]) : 



1 r 



qmit') + mm 



4 



^2g-ia;o(i+i') _^ ^t2gi<^o(t+t') 



+ {dd^ + d^d)cosujo{t-t')\, (B.IO) 
in which formula the symmetric ordering is evident. 
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